Spin Superfluidity versus Solidity of Ultracold Polar Molecules 
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We present a technique for engineering quantum magnets via ultracold polar molecules in optical 
lattices and explore exotic interplay between its spin superfluidity and solidity. The molecular 
ground and first excited rotational states are resonantly coupled by a linearly polarized microwave 
field. The spin-up (spin-down) states are presented by molecular rotational states of electric dipole 
moment along (against) the coupling field. By controlling the angle between the lattice direction and 
the coupling field, the inter-site interaction can be tuned from antiferromagnetic to ferromagnetic. 
Furthermore, the dipole-dipole interaction induces an exotic interplay between spin superfiuidity 
and solidity, and spin supersolid phases may appear in mediate regions. 

PACS numbers: 03.75.Hh, 67.85.Hj, 33.80.-b, 34.50.Rk 



In recent years, ultracold polar molecules in optical 
lattices have stimulated great interests in both funda- 
mental research and technology application (for a recent 
review, see [l| and references therein). Now, ultracold 
polar molecules such as LiCs and KRb can be prepared 
into their rovibrational ground state SiBB, 6]. Com- 
pared to ultracold atoms, ultracold polar molecules have 
some additional promising features such as their complex 
rotational and vibrational structures and strongly long- 
range anisotropic dipole-dipole interactions. The high 
controllability and tunability of internal rotational struc- 
ture of polar molecules via external DC and AC electric 
fields makes them ideal tools to design quantum logic 
gates for information processing 0, S, Sj- The strong and 
tunable dipole-dipole interaction between ultracold polar 
molecules is very different from the short-range interac- 
tions such as s-wave scattering. So that systems of ul- 
tracold polar molecules provide an excellent opportunity 
for simulating and exploring quantum phases in stro ngly 
correlated systems of long-range interactions l^, 11 , 12 |. 

The superfluid - Mott insulator phase transition, which 
is well described by a Bose-Hubbard model, has been 
experimentally demonstrated by manipulating ultracold 
atoms in optical lattices jl^. It is suggested that a super- 
soHd phase with coexistence of diagonal and off-diagonal 
long-range orders would appear in systems of off-site in- 
To reaHze the off-site interactions. 



teractions 
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15lj . To realize the off-site interactions, it 
has proposed to use dipole-dipole interactions between 
dipolar bosons [l3| or promote bosons into higher bands 
[1711 . It seems that supersolid phase could be simulated 
by spinless dipolar atoms in optical lattices, but the 
weak off-site interaction between atoms hinders further 
progress. As bosons can be represented by spins, super- 
solid states can also appear in quantum magnets, several 
numerical studies on spin-dimer models have given posi- 
tive conclusions [l8, 19, 20]. Different fro m the quantum 
magnets of conventional condensed matters, the quan- 
tum magnets of ultracold polar molecules 21 , 2^ can be 



easily tuned in experiments. 

In this Letter, we propose an alternative scheme to en- 
gineer quantum magnets of long-range dipole-dipole in- 
teractions with ultracold polar molecules in deep optical 
lattices. Usually, rotational excitations of the molecules 
are much smaller than their electronic and vibrational 
excitations. In the presence of a resonant AC electric 
field for first two rotational states, the electronic and 
vibrational excitations are almost absent, and so that 
molecular rotational states of dipole along (against) the 
electric field can be used to present the spin-up (spin- 
down) states. This spin presentation scheme differs from 
the previous ones via electronic states [2l| or first ex- 
cited rotational state in the absence of an external elec- 
tric field [i^l. Furthermore, our scheme provides a direct 
realization of spinor qubits Q rather than an equivalent 
mapping or a quantum projection. The long-range in- 
teraction between polar molecules naturally brings the 
exchange interaction between effective spins. Various 
quantum phases including superfiuid and supersolid in 
quantum spin-^ chain of one polar molecule per site are 
explored. By utilizing existed experimental techniques, 
the possibility of preparing and manipulating such kind 
of quantum magnets are briefiy discussed. 

We consider an ensemble of ultracold polar molecules 
confined in a deep optical lattice with one molecule per 
site, see Fig. [DJa). A single molecule with mass m in 
such an optical potential obeys Hamiltonian TJ^iol = 
p2/2TO -I- ijopt + Here, p'^/2m and i/°Pt are ex- 

ternal kinetic and potential energies for center-of-mass 



motions of the molecule, respectively. H 



rot _ 



the internal rotational motion characterized by the rota- 
tional constant B for electronic- vibrational ground state 
(for KRb, B « 1.1 GHz ^) and the total angular mo- 
mentum operator j. Introducing quantum numbers j and 
M for the total angular momentum and its z-component, 
the rotational eigenstates can be expressed as \j,M). 
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where erf is the z-component Pauh-matrix of molecule i, 
Tij = rijBij is the position vector between two molecules, 
and 9 is the angle between the coupling field direction 
and the lattice direction e^-, see Fig. [Hb). 

Under a Born-Oppenheimer approximation, the sys- 
tem of ultracold polar molecules in deep optical lattices 
can be described by 

-ff = ^ ^ Ucrcrriicr [Uia - 1) + U^l ^ Tli^n^l 

+ E^^^-* E(4^^- + '^-^-)-'^E^^ (3) 
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with spin indices a-{—\ or J,), lattice indices i and 



Figure 1: (Color online) (a) The schematic diagram for di- 
atomic polar molecules in an ID optical lattice, (b) An AC 
AC 

electric field E along z-axis is applied, where 9 is the angle 
between the z-axis and the lattice direction, (c) The AC field 
resonantly couples two molecular rotational states |0,0) and 
11,0), i.e., Lu = 2B. 



In the presence of an AC electric field, E^^(i) 



e^+cc, the single-molecular Hamiltonian will 



have a part H^^{t) = -d • e^E'^^e"*"* + h.c. from 
the electric dipole interaction. Due to (0,0| dz |1,0) — 
(1,0|4 |0,0) = d/V3 (in which = ■ d) [l3|, the 
AC electric field will mix |0,0) and |1,0). The frequency 
ui is chosen as 2B which is the transition frequency for 
|0,0) |1,0), see Fig.[ll[c). The AC electric field of such 
a frequency can be generated by a microwave cavity. 

To keep M unchanged, one has to use a microwave 
of linear polarization. Due to transitions |0,0) \j,0) 
of j > 2 are far from resonant, the excitations to states 
|j, 0) of J > 2 can be ignored. Thus, the molecular states 
are coherent superposition of |0,0) and |1,0). Under a 
rotating wave approximation, the Hamiltonian in basis 



{|0,0),|1,0)} reads [ll| 



^rot ^ ^AC 



n 
n 



(1) 



with n = E^^ {l,0\d;,\0,0) = E^^d/V^. Introduc- 
ing a new basis of 'dipole-up' |t) = (|0, 0) -I- |1, 0)) 

and 'dipole-down' ||) — :^(|OiO) ~ ll^O))) one can ob- 
tain H^°^ + H^^ = -\^(J'' with denoting the z- 
component PauH matrix. Obviously, (t|d|t) = dez/\/3, 
(i|d|i) = -de,/V3 and (T|d|i) = a|d|T) = 0. Thus 
the dipole-dipole interaction if^-^ between molecules i 
and j reads 



di ■ dj - 3(di • etj){dj ■ e,, 



,1-3 cos2 e 
3|r,,|3 



-cr,- a A 



(2) 



j, molecular creation (destruction) operators cj^. (cjcr), 
on-site interaction C/ctct', nearest-neighbor tunneling t, 
molecular number operators nia = c\^Cia, and z- 
component spin operator Sf — (c|^|Ci| — c|^|Cij^)/2. 

Here, H]^^ are off-site dipole-dipole interaction between 
molecules, and the effective magnetic field strength h is 
just the Rabi frequency f2. 

In the hard-core limit (t ^ Uaa' ) of one molecular per 
site {n = ni + = 1), the 'dipole-dipole' interaction is 



^dd — 



with = 4d^(l - 3 008^6') /a I (j - i)a\ 



where a is the distance between two neighboring sites. 
By adjusting 9, Jij can be tuned from positive to nega- 
tive continuously, that is, the spin-spin interaction can be 
changed from antiferromagnetic (AF) to ferromagnetic 
(FM). Due to the dipole-dipole interaction fast decays 
according to l/r^ , we only take into account the nearest- 
neighbor and the next-nearest-neighbor interactions with 
J = 4^2(1 - 3cos^0)/3a^ and J' = J/8, respectively. 
Thus the Hamiltonian |[3t is equivalent to an anisotropic 
XXZ Heisenberg model [23| 

+J'J2s^S^^-hJ2s^, (4) 

((ife)) i 

with = J + At'^/U^i - At^/Uii - 4iVt^TT and J^y = 
4iVC^Ti) the '-f' and '-' signs before Jxy correspond 
to the cases of fermionic and bosonic molecules, respec- 
tively. Here the At'^/Uaa' in Jz and Jxy are induced by 
spin-dependent collisions in second-order tunnehng pro- 
cesses. Below we shall only consider the bosonic case, 
that is, J = Jz — Jxy. The fermionic case can be mapped 
onto the bosonic case with a rotation of tt along z axis 
on one sublattice. The above XXZ Heisenberg chain (|4]) 
is equivalent to a hard-core boson model [H, [l^l 

{ij) (ij) 

+J' '^UiUk - ^J,^ni, (5) 
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Figure 2: Various magnetic states of an ID XXZ spin Hamil- 
tonian ([4]). (a) Antiferromagnetic state, (b) Intermediate 
state, (c) Spin-flop state, (d) Paramagnetic state. 



S ^ 



and /i 



with ai 
h + .], + .]'. 

The zero-temperature ground state of XXZ Heisen- 
berg Hamiltonian can be antiferromagnetic (AF), spin- 
flop (SF) or paramagnetic (PM) state ji^l, which corre- 
spond to mass density wave, superfiuid and Mott insula- 
tor states in the hard-core boson system, respectively. It 
was also suggested that a first-order phase transition be- 
tween the antiferromagnetic and spin-flop states is split- 
ted into two second-order phase transitions by an 'in- 
termediate' (IM) state, in which diagonal (S'^) and off- 
diagonal {S^y) long-range orders coexist, which corre- 
spond to the coexistence of solid and superfluid order in 
the hard-core boson system, thus it was also termed as 
supersolid. These magnetic states can be characterized 
with various sublattice magnetization as in Ref. [H], see 
Fig. [2 

The existence of diagonal long-range order is char- 
acterized by the linear increase of static structure fac- 
tor S(Q) or S(t:) with the lattice size N [IBl, where 
'S'(fc) = (I Y^i exp{ikl)S[\'^) /N is the Fourier transforma- 
tion of density-density or spin-spin correlations. While 
the existence of off-diagonal long-range order is charac- 
terized by the non-zero superfluid density or the spin- 
stiffness ps, which is deflned as the second derivative of 
the ground-state energy per site Eq with respect to a 
twist $ in the boundary condition, see Ref. MM. 

To give a qualitative picture, we use a Gutzwiller 
variational ansatz as the ground state wave function 
l^*) = riil'/'j)- The wave functions for each site 
are ex pres sed in the basis of Fock states for hard-core 
bosns [23, m, i.e., |</)^) = /^(0)|0), -t- Using 
the mean-fleld approximation, the system is divided into 
two intersecting sublattices A and B with /i(0) = sin a, 
— cosa for i & A, and fj{0) = sin/3, = cos/? 

for j G B. By minimizing the energy per site Eq — 
{"i'lHl'i/) /N with respect to a and (3, the ground states 
in J-h plane are shown in Fig. [3l In region of J > 0, 




Figure 3: The h-J phase diagram. The order parameters for 
a scan along the dashed line of J/ Jxy ~ 1 are shown in Fig.Sl 
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Figure 4: The mean-field results of static structure factors 
S{0), S{tt) and superfiuid density ps versus the effective mag- 
netic field h at J/ Jxy = 1 and J' = J/8. 



when h gradually increases. Because of the dipole-dipole 
interaction, values of J' are nonzero and supersolid states 
may appear in a certain region which increases with J. 

For a given mean-fleld ground state, one can obtain the 
static structure factors S{0)/N = (cos 2a -I- cos 2/3)^/16, 
S{tt)/N — (cos 2a — cos 2/3)^/16, and the superfluid den- 
sity Ps — \ sin 2a sin 2/3 by using similar approaches in 
[29||. In Fig. |4]we show how S{0), 8(71) and ps depend 
on h for J/ Jxy = 1 and J' = J/8. From the flgure we can 
see that in the spin supersolid or intermediate state (re- 
gion marked 'IM' in the flgure), the diagonal long-range 
order ^(Tr) and off-diagonal long-range order ps are all 
non-zero. 

We have also carried out stochastic series expan- 
sion (SSE) quantum Monte Carlo (QMC) simulation 



M[3ii,E 



. In our calculation, the inverse tempera- 
ture is chosen as 4:N. The superfluid density ps is mea- 
the phase transition AF-IM-SF-PM occurs sequentially sured via the fluctuations of the winding number. For a 
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Figure 5: The QMC results of S{ti) and pa correspond to the 
mean-field ones shown in Fig. [H The inset shows the finite 
size scaling: S{n)/N vs 1/N for different values of h. 

system of J/Jxy — 1 and J' / Jxy — 0.125, the dependence 
of S{'k) and ps on h are shown in Fig. [5l It shows that 
the QMC result is quantitatively different from the mean 
field one. First, the region of AF state is greatly reduced 
because of quantum fluctuation. Second, there is a sud- 
den drop of S{-k)/N as h increases, which shows a first 
order phase transition. We find this sudden drop even 
appears when J' j Jxy = 0.25 and 0.5. After the transi- 
tion, the scaHng behavior (see the inset of Fig. [5]) shows 
S{tt)/N is still finite, it will gradually decrease to as ft- 
increases further, thus a spin supersoHd state appears. 

Below we briefiy discuss the experimental possibility. 
It has demonstrated how to prepare polar molecules in 
optical lattices with one molecule per site and cool them 
to their rovibrational ground states [3, 0, 0, S, S|- To 
ensure sufficiently large dipole-dipole interaction, one can 
choose the heteronuclear diatom molecules such as LiCs 
and KRb which have very large permanent electric dipole 
moments in an order of a few Debye (for KRb in its 
singlet rovibrational ground state, d ~ 0.566 Debye 0). 
The resonant coupling between states |0, 0) and |1, 0) can 
be achieved by a linearly polarized microwave. The sign 
and strength of inter-site interaction can be adjusted by 
controlling the wavelength of the optical lattice and the 
angle between the lattice direction and the coupling field. 
The effective magnetic field strength can be tuned by 
varying the coupHng strength, i.e., the Rabi frequency. 

In summary, by using ultracold polar molecules in op- 
tical lattices and coupling their first two rotational states 
via an AC electric field (microwave) , we show an effective 
scheme for implementing long-range dipole-dipole inter- 
action. By controlling the optical lattice and the cou- 
pling field, a quantum spin chain of long-range exchange 



interaction can thus be simulated. Quantum phase tran- 
sitions induced by tuning the external field strength are 
explored. Experimental realization of our scheme may 
give new insights into some open questions in condensed 
matter physics including the competition between solid 
and superfiuid orders and the existence of supersolid. 
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